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^ . Abstract 



The parton content of real (P 2 = 0) and virtual (P 2 7^ 0) transverse photons 
7(P 2 ) is expressed in terms of perturbative pointlike and nonperturbative hadronic 
(VMD) components, employing recently updated parton distributions of pions and 
protons. The resulting parameter-free and perturbatively stable LO and NLO 
parton densities p^ p \x,Q 2 ) are smooth in P 2 and apply to all P 2 > when- 
ever 7(-P 2 ) is probed at scales Q 2 3> P 2 where transverse photons also dominate 
physically relevant cross sections. Predictions are given for the structure function 
F% \x,Q 2 ) and / 7 ( p2 )(x, Q 2 ), and are compared with all relevant data for real 
photons as well as with recent data for virtual photons as extracted from DIS ep 
dijet events. Simple analytic parametrizations of our predicted parton distributions 
are presented for the real photon in LO and NLO, and for the virtual photon in LO 
which, within sufficient accuracy, may be also used in NLO-QCD. 



1 Introduction 



Modern theoretical QCD studies J^, 0, [!J of the parton distributions of real , i.e. on-shell, 
photons / 7 (x,Q 2 ), / = q, q, g, agree surprisingly well with measurements of the (anti) 
quark and gluon contents of the resolved real photons as obtained from e + e~ and ep 
reactions at collider energies (for recent reviews, see ||, [5|, For clarity let us denote 
the resolved real target photon with virtuality P 2 = — p 2 ~ by 7 = 7(P 2 ~ 0) which is 
probed by the virtual probe photon 7*(Q 2 ), Q 2 = —q 2 , via the subprocess 7*(Q 2 )7 — > X as 
in e + e~ — > e + e~X. Here, p denotes the four-momentum of the photon emitted from, say, 
an electron in an e + e~ or ep collider. In the latter case it is common to use Q 2 instead of 
P 2 for denoting the photon's virtuality, but we prefer P 2 for the subprocess 7(P 2 )p — ■> X 
according to the original notation used in e + e~ annihilations. (Thus the factorization 
scale in f J ^ p2 \x,Q 2 ) refers now to some properly chosen scale of the produced hadronic 
system X, e.g. Q ~ pj, in high-p T jet events, etc.). 

In general one expects [7-12] also a virtual photon 7(P 2 7^ 0) to possess a parton 
content P (p2) {x, Q 2 ). It is a major problem to formulate a consistent set of boundary 
conditions which allow for a calculation of f 1<yP2 \x, Q 2 ) also in the next-to-leading order 
(NLO) of QCD as well as for a smooth transition to P 2 = 0, i.e. to the parton distributions 
of a real photon (see refs. [1^, [11], [L3|] for a detailed discussion). Indeed, experimental 
studies of the transition of the deep inelastic (di-)jet cross section from the real photon 
to the virtual photon region at HERA point to the existence of a nonvanishing, though 
suppressed, parton content for virtual photons |14[]. These measurements have triggered 
various analyses of the dependence of the ep jet production cross section on the virtuality 
of the exchanged photon [|l|, |16[ and experimental tests of such predictions will elucidate 



the so far unanswered question as to when a deep inelastic scattering (DIS) ep process 
is eventually dominated by the usual 'direct' 7* = ■y(P 2 ) induced cross sections, not 
contaminated by the so far poorly known resolved virtual photon contributions. More 
recently, NLO calculations of the (di-)jet rate in ep (and e^y) scattering, which properly 
include the contributions of resolved virtual photons, have become available 0] and the 



resolved virtual photon contributions have already been included in a Monte Carlo event 
generator [18| as well. 



It is the main purpose of the present paper to formulate a consistent set of boundary 
conditions, utilizing our valence-like input parton distributions at the universal target- 
mass independent [Q, [TP], ^0] low resolution scale Ql = fi 2 ~ 0.3 GeV 2 , which allow for 
a calculation of f^ p2 \x, Q 2 ) also in NLO-QCD as well as for a smooth transition to 
the parton distributions of a real photon, P 2 = 0. We shall furthermore employ the 
recently updated parton distributions of the pion |[21|| , f n (x,Q 2 ), which are required for 
describing, via vector meson dominance (VMD), the hadronic components of the photon. 
It should be noted that the pionic gluon and sea densities, g' K {x, Q 2 ) and q w (x,Q 2 ), can be 
uniquely derived |2TJ, ^2] from the experimentally rather well known pionic valence density 
v n (x,Q 2 ) and the (also recently updated |2^| dynamical) parton distributions f(x,Q 2 ) of 
the proton. Thus we arrive at essentially parameter- free predictions for f l( - p2 \x,Q 2 ) 
which are furthermore in good agreement with all present measurements of the structure 
function of real photons, F^x.Q 2 ). 

In Sec. 2 we discuss the basic theoretical framework necessary for the presentation 
of our model for the parton distributions and structure functions of real photons, and 
compare the resulting predictions with recent experiments. Sec. 3 contains the formu- 
lation of our model for the parton distributions of virtual photons, together with some 
quantitative predictions for structure functions as well as a comparison with very re- 
cent data extracted from DIS dijet events. Our conclusions are drawn in Sec. 4. In the 
Appendix we present simple analytic parametrizations of our LO and NLO predictions 
for the parton distributions of real and virtual photons. 
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2 The Parton Content of Real Photons 



High energy photons are mainly produced by the bremsstrahlung process e(k) — > e(k') 
+ 7(p). Here p 2 = (k' — k) 2 = 2m 2 — 2k ■ k! determines the virtuality of the produced 
photon which is declared as real (virtual) whenever P 2 = —p 2 is smaller (larger) than some 
Pq arbitrarily fixed experimentally (typically, P 2 ~ 10~ 2 GeV 2 ). Accordingly, the theo- 
retical analysis is usually subdivided into two distinct parts corresponding to P 2 > P 2 . 
In terms of the photonic parton distributions, the structure function F% of a real photon 



7 



l(P 2 £ Pq) is given in NLO(MS) QCD by 
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[C q ®{q + qr + 2C g ®g^} + ^e 2 C 7 (x)} (1) 



where ® denotes the usual convolution integral. Here g 7 (x, Q 2 ) = g 7 (x, Q 2 ) and g^(x, Q 2 ) 
provide the so-called 'resolved' contributions of 7 to F2, while C 1 provides the 'direct' 
contribution as calculated according to the pointlike 'box' diagram 7* (Q 2 )j(P 2 ~ 0) — > qq. 



The Wilson coefficients p4| C w are given by 
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where the convolution with the [ ] + distribution can be easily calculated using, for example, 
eq. (A. 21) of the first article in ref. [113, while the direct term in (1) is [BH 
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The appropriate coefficient function C Qj l and C g ^ for the longitudinal structure function 

'7 
L 



F7 — F% — 2xFi may be found in [TC 



In order to avoid the usual instabilities encountered in NLO(MS) in the large-x region 
due to the ln(l — x) term in (3), we follow ref. Jl[] and absorb such terms into the photonic 
MS quark distributions in (1): this results in the so-called DIS 7 factorization scheme 



which originally has been introduced for real photons by absorbing the entire C 7 term 
appearing in (1) into the NLO(MS) quark densities Q 7 (x,Q 2 ), i.e. 

with C 7 (x) given by eq. (3). How much of the 'finite' terms in (3) is absorbed into the 
MS distributions in (4), is of course arbitrary and a matter of convention ||. Since 
such different conventions || f| turn out to be of minor importance for our quantitative 
results to be discussed below, we prefer to stick to the original DIS 7 scheme [|I|] as defined 
in eq. (4). Furthermore, the redefinition of the parton densities in (4) imply that the 
NLO(MS) splitting functions kjv(x) of the photon into quarks and gluons, appearing 
in the inhomogeneous NLO RG evolution equations [|l], |Kj for p(x,Q 2 ), have to be 



transformed according to [ffl, 26 



|dIS 7 — kq I MS — e q ^qq ® ^7 

^Idis, = kV\ m -2j2e 2 q P^®C,. (5) 

q 

The LO expression for F£ is obviously entailed in eq. (1) by simply dropping all NLO 
terms proportional to C q<g as well as C 7 . 

In NLO the expression for in the above DIS 7 factorization scheme is given by 
retaining the C Qt g terms while dropping the destabilizing C 7 terms in eq. (1), which has 
already been absorbed into the quark densities according to eq. (4). Furthermore the 
nonperturbative hadronic VMD input / 7 (x, Qq) in NLO refers to the partons in the DIS 7 
scheme which guarantees the perturbative stability of the resulting F^(x, Q 2 ) provided this 
input is given by the NLO f*(x, Qq) of ||21|| , while the corresponding input in LO is given 



via VMD by the LO f n (x,Qo) of pl[ . We shall assume that the input resolution scale 
Qq = fi 2 ~ 0.3 GeV 2 for the valence-like parton structure is universal, i.e. independent of 
the mass of the considered targets p, tt, 7, etc. [19], p2| . The hadronic VMD ansatz 
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for / 7 (x, fi 2 ) is based on a coherent superposition of vector mesons |TT 



6 6 

|7) M 2 ,had ^ 7- Ip)m 2 + 7 k)^ 2 ( 6 ) 

where the 0-meson contribution is considered to be strongly suppressed at /i 2 <C mi. 
Assuming, within the operator product expansion (OPE), 

(p|0,|p>„ 2 = (^|O a |a.>^ = 27 3<? e^<p|O a |^>^ = (tt 1 Q | tt )^ (7) 

( P \o g \p)^ = Hcg^ 2 = (7t |o 3 |7t ) m2 (8) 

and (p\O g \uj) = due to isospin conservation, one obtains, 

+ M) 7 (x,/i 2 ) = a(g 2 p + gl + 2 g p g w cos 6)(u + u) n °(x, p 2 ) 

(d + dy(x,fi 2 ) = a(g 2 p + gl - 2 g p g„ cos 6) (d + d)* (x,(j, 2 ) 

(s + sy(x,^ 2 ) = a(g 2 p + g 2 J(s + sr\x,fi 2 ) = 

g\x^ 2 ) = a{g 2 p + gl)g-\x^ 2 ). (9) 

Here O q g refer to the leading twist-2 quark and gluon operators in the OPE formalism, 
9' = (7|0 9 |7)had etc., and g v = 4tt 7 ' f v with 

£ = 0.50, gl = 0.043, (10) 

i.e. f 2 /4:7T = 2.0 and / 2 /47r = 23.26, as obtained from a zero-width calculation of the 



relevant leptonic widths r(V" — > = a 2 mvgv/3 presented in p7|. The omission of a 

finite-width correction for g 2 is due to the central role J?7|] of the precise results in [28] 
which do not require such a correction in contrast to the situation for the less precise 
resonance analysis at e + e~ colliders |2^ . 



For the a priori unknown coherence factor (fit parameter) cos# in eq. (9) we take 
cos# = 1, i.e. we favor a superposition of u and d quarks which maximally enhances 
the contributions of the up-quarks to F^ in eq. (1). This favored value for cos 9 is 
also supported by fitting cos 6 in (9) to all presently available data on F^x, Q 2 ), to be 
discussed below, which always resulted in cos# ~ 1 in LO as well as NLO. This is also in 
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agreement with the LO results obtained in ref . [II] . The LO and NLO input distributions 
f n (x,fi 2 ) of the pion in (9) are taken from a recent analysis |21] which correspond to 

13 H Vlo = °- 26 GeV2 and /4lo = °- 40 GeV2 in L0 and NL ° ; respectively. Since 
by now all free input quantities have been fixed in eq. (9), we arrive at rather unique 
parameter- free predictions for / 7 (x, Q 2 ) and F^x, Q 2 ). 

The calculation of f 1 {x,Q 2 ) at Q 2 > fi 2 follows from the well known inhomogeneous 
RG evolution equations in LO and NLO (see, for example, [1], ^(J) which we solve, as 
usual, analytically for the n-th Mellin moment of f J (x, Q 2 ), followed by a straightforward 
Mellin-inversion to Bjorken-x space. The explicit formal solutions can be found in eqs. 
(2.12) and (2.13) of the first article of ref. [ID. The general structure of these solutions is 



P(x,Q 2 ) = r pl (x,Q 2 ) + fl d (x,Q 2 ). 



11) 



Here f pl denotes the perturbative 'pointlike' solution which vanishes at Q 2 = /i 2 and is 
driven by the pointlike photon splitting functions k^'(x) appearing in the inhomoge- 
neous evolution equations, while /^ ad depends on the hadronic input / 7 (x, /i 2 ) in eq. (9) 
and evolves according to the standard homogeneous evolution equations. We treat these 
solutions in precisely the same way as discussed in detail in JTU |, except that we implement 
an improved treatment of the running a s (Q 2 ) by exactly solving in NLO (MS) 



da s {Q 2 
dlnQ 2 



A) 2r n 2\ _A>_ ^3/,o2\ 



(12) 



numerically |23] using a s (M|) = 0.114, rather than using the usual approximate NLO 



solution (e.g. eq. (14) in [|10j) which becomes sufficiently accurate only for Q 2 ^ m 2 ~ 
2 GeV 2 |1. Here, (3 = 11 - 2//3 and /3i = 102 - 38//3. 

The prescription for the VMD ansatz in eq. (9) at the input scale /i 2 , together with 
cos 9 = 1 as discussed above, yields a simple expression for the general Q 2 -dependence of 

p(x,Q 2 y. 



P(x, Q 2 ) = f pl (x, Q 2 ) + a Gj P(x, Q 2 ) + 5 f -(G 2 U - G d ) 8 *(x, Q 2 ) 



(13) 



with S u = —1, 3d = +1 and 5 S = 5„ = 0, and where the index 7r obviously refers to ir° and 
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Gl 



G ' 



0.836 
0.250 



Gl 



g 2 + gl = 0.543. 



(14) 



Simple analytic LO and NL0(DIS 7 ) parametrizations for the pointlike piece fZ(x,Q 2 
are given in the Appendix, whereas the ones for /""(x, Q 2 ) can be found in |^T] . 



The photonic quark distributions discussed thus far and which appear in eq. (1) are 
adequate for the / = 3 light u, d, s flavors. Since heavy quarks h = c, b, t will not be 
considered as 'light' partons in the photon (as in the case of the proton p3| and the 



pion pT| ; p2]]), their contributions to F 2 have to be calculated in fixed order perturbation 
theory according to the 'direct' box-diagram 7*(Q 2 )7 — > hh expression, i.e. the usual 
Bethe-Heitler cross section GO 



1 



x 



F 2,h( X 'Q 2 ) 



2>ei-9{(5 2 ) \(3 
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x 2 + (1 - x) 2 + x(l - 3x) —r± - x 



2 Sm A h 



In 



1-/3 



(15) 



Q 2 ~ Q 4 

where (3 2 = 1 — Am 2 h /W 2 = 1 — 4m^x/(l — x)Q 2 . This expression has to be added to 
eq. (1) and a similar expression holds for the longitudinal structure function [ TD | = 
F 2 — 2xFi. The 'resolved' heavy quark contribution [K| to F 2 in (1) has to be calculated 
via 7*(Q 2 )# 7 -> hh, 

- 1 dz 



z \z 



(16) 



where \ fl (Q )gJ ^ hh (x, Q 2 ) is given by eq. (15) with e\a -> e 2 h a s (/j, 2 F )/Q, z mill = 
x(l + Aml/Q 2 ) and y? F ~ 4m 2 ^l[|. This 'resolved' LO contribution should be included 
in eq. (1) as well. To ease the calculations we shall keep these LO expressions (15) and 
(16) also in NLO, since the full NLO expressions for heavy quark production |32| turn out 
to be a small correction to the already not too sizeable (at most about 20%) contribution 
in LO. Notice that such small corrections are not larger than ambiguities due to different 



choices for and for the factorization scale //p. For our purposes it is sufficient to 
include only the charm contributions which will be calculated using m c — 1.4 GeV. 



Having outlined the theoretical basis for our photonic parton distributions, we now 
turn to the quantitative results. First we apply our parameter-free predictions for 
/ 7 (x,Q 2 ) to the structure function of real photons which, according to eq. (1) and the 
above results, is finally given by 

-F^(x,Q 2 ) = 2 £ e 2 q \q\x,Q 2 ) + [C q ®q^ + C g ®g^] 

X j I 27T 

q=u,d,s \ 

+ Ii^ c (*,Q 2 ) + -F( c (x,Q 2 )} (17) 

X X J 

where p{x, Q 2 ) refers to the DIS 7 factorization scheme defined in (4) and the charm 
contributions ^ and F!f c are given by eqs. (15) and (16), respectively. In fig. 1 we 
compare our LO and NLO predictions with all available relevant data [^] for i 7 ^ of the 



real photon. Our present new NLO results are rather similar to the ones of AFG 0, 
but differ from our previous (GRV 7 ) predictions [|T| which are steeper in the small-x 
region, as shown in fig. 1, because the dominant hadronic (pionic) sea density q n (x, Q 2 ) is 
steeper since it has been generated purely dynamically from a vanishing input at Q 2 = fi 2 
0, ^Of. Similarly the SaS ID expectations, for example, fall systematically below the 



data in the small to medium Q 2 region around Q 2 ~ 5 GeV 2 , partly due to a somewhat 
different treatment of the hadronic coherent VMD input as compared to our results in 
eqs. (9) and (13) and (14). The relevant LO and NLO photonic parton densities are 
compared in fig. 2 at Q 2 = 10 GeV 2 . For illustration we also show the purely 'hadronic' 
component (homogeneous solution) in (11) of f 1 which demonstrates the dominance of 
the 'pointlike' component (inhomogeneous solution) in (11) for u 1 and d 1 in the large-x 
region, x > 0.1. In fig. 3 we show our predictions for xu l {x,Q 2 ) and xg 1 {x,Q 2 ). The 
parton distributions of the photon behave, in contrast to the ones of a hadron, very 
differently in the limits of large and small x. In the former case, the purely perturbative 
pointlike part in (11) dominates for x £ 0.1, especially for the quark distributions. On 
the other hand, this uniquely calculable contribution amounts at most to about 20% at 
very small x where the hadronic VMD component in (11) dominates, giving rise to a 
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very similar increase for x — > as observed in the proton case. In fig. 3 we also show 
our valence-like inputs at Q 2 = yU 2 Q NLO which become (vanishingly) small at x < 10 -2 . 
This illustrates the purely dynamical origin of the small-x increase at Q 2 > fi 2 . Also 
noteworthy is the perturbative LO/NLO stability of m 7 (x, Q 2 ) which is almost as good as 
the one required for a physical quantity like F^x, Q 2 ) in fig. 1. The situation is, as usual 
0, [19], |23f , different for g^^x, Q 2 ). Nevertheless, despite the sizable difference between the 
LO and NLO gluon distributions in fig. 3 in the small-x region, the directly measurable 
F% and the gluon-dominated heavy quark contribution in eq. (16) shows a remarkable 
perturbative stability p3 |. Finally, we compare in fig. 4 our predictions for xg' y (x,Q 2 ) 



at Q 2 = (pq?) 2 — 75 GeV 2 with recent HERA (HI) measurements |34j. Our somewhat 
flatter results for xg 1 in the small-x region, as compared to the older GRV 7 expectations 
0, is caused by the recently favored flatter gluon distribution in the proton |£J which 
determines g 1 via g n cf. eq. (9), at small x. 

Finally it is interesting to consider the total momenta carried by the photonic partons, 

M](Q 2 )= E t xP(x,Q 2 )dx. (18) 



Inspired by the ideas and suggestions put forward in refs. [35, 11], it has been conjectured 



recently [^F|| that this leading twist-2 quantity Ml should satisfy, in LO-QCD, 

M2{Q 2 ) ~ H,(g 2 ) (19) 

where the well known dispersion relation relates the hadronic part of the photon's vacuum 
polarization 

^ Q 2 r °h(s) 



n ^'^y 4 „ ; if#* (2o) 

to (Tft = cr(e + e~ — > hadrons). It should be noted that Uh{Q 2 ), being an experimental 
quantity, includes, besides the usual twist-2 term, all possible nonperturbative higher- 
twist contributions. The 'consistency' relation (19) is, however, expected to hold already 
at Q 2 >, 2 to 4 GeV 2 to within, say, 20 to 30% where the twist-2 component in Tl^iQ 2 ) 
may become dominant, as possibly indicated by DIS ep processes. Indeed, our LO results 
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imply M 2 7 (2 GeV 2 )/a ~ 0.976 and M 2 7 (4 GeV 2 )/a ~ 1.123 which compares favorably 
with |7j Il h (2 GeV 2 )/a = 0.694 ±0.028 and n fc (4 GeV 2 )/a = 0.894 ±0.036, respectively. 



3 The Parton Content of Virtual Photons 

Next we turn to the somewhat more speculative concept and models of 'resolved' virtual 
photons (P 2 7^ 0). The real photons considered in Section 2 are those whose virtuality P 2 
is very small, i.e. of the order P^ in = 0(m 2 ) or, experimentally, at least P 2 < P 2 ~ 10~ 2 
GeV 2 . The flux of virtual photons produced by the bremsstrahlung process e(k) — > 
e(k') + 7(p), P 2 = -p 2 = ~(k' - k) 2 > P 2 is given by @ 

'l + il-yf 1 2m 2 y 



a 



y P 2 P 4 



(21) 



L _ a 2(1 -y) 1 

7 7 (pa)/eW - ^ y ~pi ^ ll > 

with y = E y /E e and T(L) denoting transverse (longitudinal) photons. Whenever these 
virtual photons, with their virtuality being entirely taken care of by the flux factors in 
(21) and (22), are probed at a scale Q 2 ^> P 2 they may be considered as real photons 
which means that 



(i) effects due to fb P 2y e should be neglected since the corresponding longitudinal cross 
sections are suppressed by powers of P 2 /Q 2 ; 

(ii) cross sections of partonic subprocesses involving 7(-P 2 ) should be calculated as if 
P 2 = due (partly) to the P 2 /Q 2 power suppressions of any additional terms. 



This latter rule implies in particular that the NLO 'direct' contribution C y /p2\(x) to 
F^ P \x,Q 2 ) has to be the same C 7 (x) as for real photons in eqs. (1) and (3), i.e. has 
to be inferred from the real photon subprocess 7*(Q 2 )7 — > qq, and not from the doubly- 
virtual box 7*(Q 2 )7(-P 2 ) — > qq as originally proposed 0, § and used |10| . Thus we can 



implement the same DIS 7 factorization scheme as for real photons in eq. (4), and the 
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structure function of virtual photons becomes formally very similar to eq. (17): 
l -Fr 2 \x,Q 2 ) = 2 £ e 2 {g^(*,Q 2 ) + ^ [C q ®<f^ + C g ®g^ 

q=u,d,s \ 

+ -FI c (x,Q 2 )+-FC (P2 \x,Q 2 )} (23) 

X X J 

with C q ^ g {x) given in (2). The 'direct' heavy (charm) quark contribution is given by 

eq. (15) as for real photons since F^f ^ = F^ c as follows from our consistent strict 

~i( p2 ) 

adherence to point (ii) above. The 'resolved' charm contribution Pf c is as in eq. (16) 
with the gluon distribution g J (z,fip) — > g l( - p2 \z , fip) . 

The above consistency requirements afford furthermore the following boundary con- 
ditions for f l( - p2 \ cf. eq. (13), 

f^ 2 \x, Q 2 = P 2 ) = f^ 2 \x, P 2 ) = V (P 2 )fUx, P 2 ) (24) 

in LO as well as in NLO. Here P 2 = max (P 2 , jj?) as dictated by continuity in P 2 [l(| and 



fj(P 2 ) = (1 + P 2 /m 2 )~ 2 is a dipole suppression factor with m 2 = 0.59 GeV 2 . The second 
equality in eq. (24) follows from the consistency requirement C 7 (p2) = C 1 and consequently 
the application of the same MS — > DIS^ factorization scheme transformation as for the 
real photon, cf. eq. (23). The scale P 2 is dictated not only by the above mentioned 
continuity requirement, but also by the fact that the hadronic component of f l( ' p2 \x^ Q 2 ) 
is probed at the scale Q 2 = P 2 || [H], |Tl| where the pointlike component vanishes by 
definition. The boundary condition in eq. (24) guarantees, as should be evident, a far 
better perturbative stability as compared to the situation in [|K| where the NLO input 
differed drastically from its LO counterpart (cf. eq. (8) in ref. [|l0| ). 

The evolution to Q 2 > P 2 is now analogous to the case of real photons in the previous 
section and the general solution for the resulting parton distributions is similar to the one 
in eq. (11) and eq. (13), 

f^ p2 \x,Q 2 ) = fZ ip2) ^Q 2 ) + fS\x,Q 2 ) (25) 



M p2 )i 



ii 



G 2 f r(x,Q 2 ) + 5 f l -{G 2 u -G 2 d ) s-(x,Q 2 



with 5f as in eq. (13) and where f^ P refers again to the solution of the homogeneous RG 
evolution equations, being driven by the hadronic input in (24), which is explicitly given 
by eq. (13) of ref. [|Kj. Its parametrization is fixed by the available parametrization [^l[ 
for f n (x,Q 2 ) in (25). The inhomogeneous 'pointlike' solution in (25) is explicitly given 
by eq. (12) of ]nj where L = a s (Q 2 )/a s (P 2 ). A parametrization of f$ p \x, Q 2 ) in LO is 
thus easily obtained from the one for the real photon fZ(x, Q 2 ) in (13) in terms of In L^ 1 = 
In [a s (fi 2 )/a s (Q 2 )], where now a s (/i 2 ) has simply to be replaced by a s (P 2 ) as described 
in detail in Appendix 1. Furthermore, since our NLO predictions for / 7 ^ p2 ^(x, Q 2 ) turn 
out to be rather similar to the LO ones, as will be shown below, the simple analytic 
LO parametrizations for p( p2 \x, Q 2 ) can be used for NLO calculations as well. This is 
certainly sufficiently accurate and reliable in view of additional model ambiguities inherent 
in the parton distributions of virtual photons. 

It should be emphasized that the RG resummed results in (25) are relevant when- 
ever P 2 <C Q 2 , typically [|H], [L5| P 2 ~ Q 2 , so as to suppress power-like (possi- 
bly higher twist) terms (P 2 /Q 2 ) n which would spoil the dominance of the resummed 
logarithmic contributions and, furthermore, to guarantee the dominance of the trans- 
verse photon contributions (21) to physical cross sections. For P 2 approaching Q 2 , the 
e + e _ — > e + e~X reaction, for example, should be simply described by the full fixed order 
box 7*(Q 2 )t(-P 2 ) — > qq keeping all (P 2 /Q 2 ) n terms. Since the full perturbative 0(a s ) 
corrections to this virtual box have not been calculated yet, it is not possible, for the time 
being, to determine reliably at what values of P 2 (and possibly x) this 0(a s ) corrected 
virtual box becomes the more appropriate and correct description. Similar remarks hold 
for a DIS process ep — > eX, i.e. 7(P 2 )p — > X, where 0(a s ) corrections to pointlike vir- 
tual 7(P 2 )-parton subprocesses have to be analyzed in detail in order to decide at what 
P 2 these pointlike expressions become the more appropriate description and the virtual 
photonic parton distributions (i.e., resummations) become irrelevant. 

Our strict adherence to the above point (ii) implies that the 'direct' photon con- 
tribution to any process whatsoever should always be calculated as if this photon is 
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real apart from the fact that its flux should be evaluated according to eq. (21) with 
P 2 ytz o |T2| , 15| , [16 1 . This differs from the somewhat inconsistent procedure adopted by 



the HERA-HI collaboration |14j where exact eq — > eqg and eg — > eqq matrix ele- 
ments were used for the direct photon contribution to the dijet cross section. As long as 
P 2 < jqQ 2 , the exact treatment of matrix elements, however, should not differ too much 



17j from the more appropriate treatment described above. To conclude let us stress that 
the strict adherence to point (ii), as illustrated by the foregoing examples, is not a free 
option but a necessary consistency condition for introducing the concept of the resolved 
parton content of the virtual photon as an alternative to a non-resummed fixed order 
perturbative analysis at P 2 ^ 0. This consistency requirement is related to the fact that 
all the resolved contributions due to P^^ix^Q 2 ) are calculated (evoluted) as if these 
partons are massless [7 — 11] (i.e. employing photon splitting functions for real photons, 
etc.) in spite of the fact that their actual virtuality is given by P 2 ^ 0. Thus the direct 
photon contribution should obviously be also treated accordingly. 

Now we turn to our quantitative predictions and first compare our results with the old 
PLUTO measurements []3I| for -F e ^ P \x, Q 2 ) = F^ P ^ + \F]} P ■* in fig. 5. Our very similar 



LO and NLO results, which are dominated by the almost unique 'pointlike' contribution 
in (25), are in full agreement with the limited poor statistics of the PLUTO data. The 
GRS (l(J expectations turn out to be very similar to our present ones shown in fig. 5. 
For illustration the naive (i.e. not resummed) LO quark-parton model 'box' expectation 
is shown as well by the dotted curve in fig. 5 which is given by 



1 *$£\x,Cf) = 3 £ et-[x 2 + (l-x) 2 }ln^ + ±FZf\x,Q 2 ) (26) 

X , 7T 1 J P z X 

q=u,d,s 

-F]}Q{x,Q 2 ) = 3 £ et-4x(l-x) + ±F2f \x,Q 2 ) (27) 

X q=u,d,s 71 X 

with the heavy quark (charm) contribution F^f ^ given by eq. (15) and 



2m 2 1 + p 



(3x(l — x) — x 2 In 



Q 2 l-P 



:28) 
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More detailed predictions for \x,Q 2 ) are presented in figs. 6a and 6b for various 
virtualities P 2 and scales Q 2 . Since the 'pointlike' component in (25) is uniquely calculable 
perturbatively, a detailed measurement of the x and P 2 dependence at various fixed values 
of Q 2 , as shown in figs. 6a and 6b, would shed light on the theoretically more speculative 
and far less understood nonperturbative 'hadronic' contribution in eq. (25) and eventually 
establish the absolute perturbative predictions. Our LO and NLO predictions in figs. 6a 
and 6b show a remarkable perturbative stability throughout the whole x-region shown, 
except perhaps for P 2 ^> 1 GeV 2 where the perturbatively very stable [21, 22] 'hadronic' 



component in (25) becomes strongly suppressed with respect to the 'pointlike' solution 
which is less stable in the small x region, x < 10~ 2 , as is evident from fig. 6b. 

The individual LO and NLO parton distributions of the virtual photon at Q 2 = 10 
GeV 2 are shown in fig. 7a where they are compared with the ones of GRS |fL0"l . The LO 
SaS expectations [jnj are compared with our LO predictions in fig. 7b. In figs. 8 and 9 
we show our predictions for xu 1<yP2 \x,Q 2 ) and xg l( - p2 \x,Q 2 ) with particular emphasis 
on the very small x region. For comparison we also show the results for a real (P 2 = 0) 
photon. Plotting the 'hadronic' component in (25) separately in fig. 8 demonstrates that 
the perturbative 'pointlike' component in (25) dominates for x > 1CT 2 . Furthermore 
the expected perturbative stability of our present LO and NLO predictions is fulfilled. 
This is in contrast to the GRS results which are unstable |10[ throughout the whole 
x-region for P 2 >, 1 GeV 2 , as illustrated in fig. 9 at Q 2 = 100 GeV 2 , due to the very 
different perturbative (box) input in LO and NLO \IU\. In general, however, as soon 



as the perturbatively very stable 'hadronic' component in (25) becomes suppressed for 
P 2 ^> 1 GeV 2 , the remaining perturbatively less stable 'pointlike' component destabilizes 
the total results for q 1 ^ p2 \x, Q 2 ) in the very small x region, x ^ 10~ 3 , as can be seen in 
fig. 9 for w^ p2 ) at Q 2 = 100 GeV 2 (cf. fig. 6b). 

Finally in fig. 10 we confront our LO predictions for f l( - p2 \x,Q 2 ) with the effective 
parton density 

P {p2 \x,Q 2 )= £ ( q ^+q^)+ 9 - g ^ (29) 

q=u,d,s 

14 



extracted in LO from DIS dijet data by the HERA-HI collaboration (14j very recently. 
The predicted dependence on the photon's virtuality P 2 at the scale Q 2 = (pt*) = 85 
GeV 2 agrees reasonably well with the measurements in the relevant kinematic region 
P 2 <C Q 2 . This is also the case at other scales Q 2 = fp^*) and fixed values of x [|n 



not shown in fig. 10. As discussed above, it should be kept in mind, however, that for 
larger values of P 2 approaching Q 2 , which refer to the dashed curves in fig. 10, the whole 
concept of RG resummed parton distributions of virtual photons is not appropriate any- 
more. Since the resolved contributions of a virtual photon with virtuality as large as 
P 2 = 10 — 15 GeV 2 are by a factor of about 10 smaller than the ones of a real (P 2 = 0) 
photon, it is reasonable to conclude from fig. 10 that for P 2 >, 10 GeV 2 the DIS ep —>■ eX 
process considered is dominated by the usual direct 7* = 7(-P 2 ) exchange cross sections 
and not 'contaminated' anymore by resolved contributions. This furthermore explains the 



trend of the discrepancies between the data and our as well as other [11], [L2| predictions 
which can be traced to the fact that the direct photon contributions were not calculated as 
if 7*(P 2 ) was real, as required by our consistency condition (ii). Thus the direct photon 
contribution was likely underestimated, particularly at the larger value of x, x — 0.6, 
resulting in an overestimate of xf l( - p2 \x,Q 2 ) at P 2 >, 5 GeV 2 . 



4 Summary and Conclusions 



The main purpose of the present paper was to formulate a consistent set of boundary 
conditions which allow for a perturbatively stable LO and NLO calculation of the photonic 
parton distributions / 7 ^ p2 - ) (x, Q 2 ) as well as for a smooth transition to the parton densities 
of a real (P 2 = 0) photon. Employing the recently updated pionic distributions 
f' K {x 1 Q 2 ) 1 required for describing, via VMD, the nonpointlike hadronic components of 
a photon, we arrive at essentially parameter- free predictions for f^ p2 \x,Q 2 ) which are 
furthermore in good agreement with all present measurements of the structure function 
F^x^Q 2 ) of real photons 7 = 7(P 2 = 0). It should be noted that the experimentally 
almost unconstrained pionic gluon and sea distributions, g n (x,Q 2 ) and g 7r (x,(5 2 ), have 
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been uniquely derived E2| from the experimentally rather well known pionic valence 
density v n (x, Q 2 ) and the (also recently updated |23| dynamical) parton distributions of 
the proton. We have furthermore implemented these hadronic components by using a 
VMD ansatz for a coherent superposition of vector mesons which maximally enhances 
the contributions of the up-quarks to as favored by all present data. Since these 
hadronic contributions are generated from the valence-like input parton distributions at 
the universal target-mass independent low resolution scale Qq = jj 2 ~ 0.3 GeV 2 , we 
arrive, at least for real (P 2 = 0) photons, at unique small- x predictions for x £ 10 -2 at 
Q 2 > fi 2 which are of purely dynamical origin, as in the case of hadrons. Furthermore, 
since our universal input scale /i 2 fixes also uniquely the perturbative pointlike part of 
the photonic parton distributions, which dominates for x £ 0.1, the large-x behavior of 
photonic structure functions is unambiguously predicted as well. 

Our expectations for the parton content of virtual (P 2 ^ 0) photons are clearly more 
speculative, depending on how one models the hadronic component (input) of a vir- 
tual photon. The latter is usually assumed to be similar to the VMD input for a real 
photon, times a dipole suppression factor which derives from an effective vector-meson 
P 2 -propagator, cf. eq. (24). Whenever a virtual photon is probed at a scale Q 2 ^> P 2 , 
with its virtuality being entirely taken care of by the (transverse) equivalent photon flux 
factor, it has to be considered as a real photon in the sense that cross sections of subpro- 
cesses involving 7(P 2 ) should be calculated as if P 2 = 0. In other words, the treatment 
and expressions for f y ^ p2 \x,Q 2 ) as on-shell transverse partons obeying the usual RG 
Q 2 -evolution equations (with the usual splitting functions of real photons, etc.) dic- 
tate an identification of the relevant resolved sub-cross- sections p( p2 ^X -> X' with that 
of the real photon, a(f^ p2 ^X — > X') = a(f J X — > X'). In particular, the calculation 
of F2 {p2 \x,Q 2 ) requires the same photonic Wilson coefficient C 7 (x) as for P 2 = 0, in 
contrast to what has been originally proposed |7], [J. This allows to formulate similar 
boundary conditions in LO and NLO which give rise to perturbatively stable parton dis- 
tributions, cross sections (i.e. also structure functions) of virtual photons 7(P 2 ) as long 
as they are probed at scales Q 2 ^> P 2 where Q 2 = 4m 2 , (p^*) , etc., and typically 
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P 2 < jq Q 2 . It should be emphasized that only in this latter kinematic range P 2 <C Q 2 
is the whole concept of RG resummed parton distributions of (resolved) virtual photons 
appropriate and relevant. Parton distributions of virtual photons extracted recently from 
DIS ep dijet data are in good agreement with our (parameter-free) predictions. 

Finally, we present simple analytic parametrizations of our predicted LO and 
NL0(DIS 7 ) parton distributions of real photons. From these LO parametrizations one 
can easily obtain also the ones for a virtual photon which, whithin sufficient accuracy, 
may also be used in NLO. Our NL0(DIS 7 ) parametrizations of the parton densities of the 
real photon can be easily transformed to the MS scheme according to eq. (4) which might 
be relevant for future NLO analyses of resolved photon contributions to hard processes 
where most NLO subprocesses have so far been calculated in the MS scheme. 

A FORTRAN package containing our most recent parametrizations can be obtained 
by electronic mail on request. 
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Appendix 



Simple analytic parametrizations in LO and NLO of the ' hadronic' piece of the real and 
virtual photonic parton distributions, being proportional to f n (x, Q 2 ) in eqs. (13) and (25), 
respectively, are already known according to the recently published ]2l| parametrizations 
for f n (x, Q 2 ). Therefore we only need to parametrize the remaining 'pointlike' components 
in eqs. (13) and (25). 



1. Parametrization of LO 'pointlike' photonic parton distributions 

In LO the Q 2 dependence of the 'pointlike' fp t (x,Q 2 ) term in (13) enters, apart from an 
overall l/a s (Q 2 ) factor, merely via the combination L = a s (Q 2 ) / a s (//| ) as is evident, for 
example, from eq. (2.12) of the first article in ref. []!]. Therefore, we prefer to parametrize 
the quantity fpi(x, Q 2 ) in terms of 

In [Q 2 /(0.204 GeV )2] 



s = In 



in[^ o /(0.204 GeVy 



(A.l) 



where fi^Q = 0.26 GeV 2 , which will later provide us a parametrization also for the 
virtual 'pointlike' component in (25). Our resulting 'pointlike' distributions in eq. (13) 
can be expressed by the following simple parametrizations, valid for 0.5 Q 2 £ 10 5 GeV 2 
(i.e. 0.31 < s < 2.2) and 10" 5 < x < 1 : 



a 



Q 2 



If) — 

4vr (0.204 GeV) 2 



s a x a {A + B^+Cx b 



+s Q 'exp -E + \ E'sPln- 



x 



1 — X 



where for 



a 

a 
A 
C 
E 



u 



pi 



pi - 

0.897, 
0.413, 

2.137-0.310 v^, 
-0.785 + 0.270 y/s, 
0.252 -0.065 y/s, 



a' 

b 

B 

D 



6.749 + 2.452 s - 0.226 s 2 , E' 



2.626, 

-1.049 + 0.113s, 

0.650 - 0.146 s, 

-0.116 + 0.403s -0.117s 2 

1.994 s -0.216 s 2 , 



(A.2) 



(A.3) 
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for fpi = dli = dli = s Pi = » 



a = 


1.084, 




a' 


= 2.811, 




a = 


0.960, 










a = 


0.914, 




b 


= 3.723- 


0.968 s 


A = 


0.081 - 


0.028 y/s, 


B 


= -0.048 




C = 


0.094- 


0.043 y/s, 


D 


= 0.059 + 


0.263 s 


E = 


6.808 + 


2.239 s -0.108 s 2 , 


E' 


= 1.225 + 


0.594 s 



0.073 s 2 



(A.4) 



and for 



i 
9pi 



a 


= 1.262, 


a' 


= 2.024, 




= 0.770, 






a 


= 0.081, 


b 


= 0.848 


A 


= 0.012 + 0.039 v/s, 


B 


= -0.056 - 0.044 s 


C 


= 0.043 + 0.031s, 


D 


= 0.925 + 0.316 s, 


E 


= 3.129 + 2.434 s -0.115 s 2 , 


E' 


= 1.364+ 1.227 s- 



0.128 s 2 



(A.5) 



With these parametrizations at hand in terms of s in (A.l), the appropriate ones for 
the 'pointlike' distributions /j/ P \x, Q 2 ) of a virtual photon appearing in eq. (25) are 
simply given by the same expressions above where in (A.l) /i 2 has to be replaced by 
P 2 = max(P 2 , /i 2 Q ). As discussed in Sect. 3, these parametrizations can, within sufficient 
accuracy, also be used for the parton distribution of virtual photons in NLO. 



2. Parametrization of NLO 'pointlike' photonic parton distributions 

In NLO the Q 2 dependence of the 'pointlike' distributions of real photons in (13), 
fp[(x,Q 2 ), can be easily described in terms of the following 'effective' logarithmic ratio 

MQ 2 /(0.299GeV) 2 ] 
In [/4lo/(0.299 GeV) 2 ] ^ 

to be evaluated for /Xn LO = 0.40 GeV 2 . Our NL0(DIS 7 ) predictions can now be 
parametrized as the LO ones and are similarly valid for 0.5 £ Q 2 £ 10 5 GeV 2 (i.e. 0.14 
£ s <, 2.38) and 10~ 5 £ x < 1. For convenience we include now the NLO a s in the r.h.s. 
of (A.2), i.e. 
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s a x a {A + B v ^+Cx b ) 



+s a 'exp (— E + J E's 13 In ^ 



{l-x) D 



where for f vl = u J pl 



u 



P i 



a = 1.051, 

(3 = 0.970, 

a = 0.412- 0.115 Vs, 

A = -0.028 + 0.019 s 2 , 

C = 6.726 - 3.264 v/s -0.166 s 2 

E = 4.122 + 3.170 s - 0.598 s 2 , 



a 



= 2.107, 



b = 4.544 - 0.563 s, 

B = 0.263 + 0.137 s, 

D = 1.145 - 0.131s 2 , 

E' = 1.615 s -0.321s 2 , 



for P 



P i 



a P l 



a P l 



'pi 



a = 1.043, 

(3 = 0.457, 

a = 0.416 -0.173 v/s", 

A = -0.010 v^ + 0.006 s 2 , 

C = 1.577- 0.916 y/s, 



a 

b 

B 

D 



E = 5.240 + 1.666 s - 0.234 s 2 , E' 
and for f£ = 



a = 0.901, 

(3 = 1.666, 

a = 0.844 - 0.820 

A = 0.194 

C = 0.330 -0.177 s, 



a 

b 

B 

D 



E = 2.895 + 1.823 s -0.441s 2 , E' = 



= 1.812, 

= 4.489 - 0.827 s, 

= 0.064 + 0.020 s 

= 1.122 - 0.093 s- 0.100 s 2 

= 1.284 s 



= 1.773, 



2.302 - 0.474 s 
-0.324 + 0.143 s 
0.778 + 0.502 s -0.154 s 2 
2.344 - 0.584 s. 
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Figure Captions 



Fig. 1 Comparison of our radiatively generated LO and NL0(DIS 7 ) predictions for 
F2(x,Q 2 ), based on the valence-like parameter-free VMD input in eq. (9), with 
the data of ref. |33| . For our comparison the GRV 7 HJ results are shown as well. In 



both cases, the charm contribution has been added, in the relevant kinematic region 
W > 2m c , according to eqs. (15) and (16). 

Fig. 2 Comparison of our predicted LO and NL0(DIS 7 ) distributions u 1 = m 7 , d 1 = d 1 
and gi at Q 2 = 10 GeV 2 with the LO/NLO GRV 7 densities 0, the LO SaS ID and 
2D [IT | and the NLO AFG P distributions. The 'hadronic' (pionic) components of 



our total LO and NLO results in eq. (11) are displayed by the dashed curves. 

Fig. 3 Detailed small-x (as well as large-a;) behavior and predictions of our radiatively 
generated u 1 = u 1 and g 1 distributions in LO and NL0(DIS 7 ) at fixed values 
of Q 2 . The dashed-dotted curves show the hadronic NLO contribution /^ ad to 
P = /pi + /had m ec L- (11)- The valence-like inputs at Q 2 = /^lo.nlo' according to 
eq. (9), are shown by the lowest curves referring to fi 2 . For comparison we show the 
steeper NLO GRV 7 |I] expectations as well. The results have been multiplied by 
the number indicated in brackets. 

Fig. 4 Comparison of our LO and NLO predictions for xg~ f at Q 2 = ((p^) 2 ) = 75 GeV 2 
with HERA(Hl) measurements [34]. The GRV 7 and SaS expectations are taken 



from refs. JT| and JTTJ, respectively. 

Fig. 5 LO and NLO predictions for F^ P "* = F^ P ^ + \F]} P \ The charm contributions 
have been calculated according to eqs. (15), (16) and (28), and the 'resolved' charm 
contribution to F]} P ^ is calculated analogously as in (16) wit\if2* iQ2)9 ^ h ~ h (x,Q 2 ) 
being given by eq. (28) with e\a e 2 a s ( / u 2 7 )/6. The 'box' is defined by eqs. (26) 
and (27). The PLUTO data are taken from ref. 
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Fig. 6a LO and NLO predictions for the x dependence of the virtual photon structure 
function ^ at Q 2 = 10 GeV 2 and various fixed values of P 2 <C Q 2 , neglecting 



any heavy quark contribution. For comparison we also show the NLO GRS 
results as well as the predictions for a real (P 2 = 0) photon. Notice that the dotted 
curve for P 2 = referred to as GRS obviously coincides with the GRV 7 result |jlj|. 
The results have been multiplied by the number indicated in brackets. 

Fig. 6b As in fig. 6a but at Q 2 = 100 GeV 2 . 

Fig. 7a Comparison of our predicted LO and NL0(DIS 7 ) distributions of the virtual photon 
at Q 2 = 10 GeV 2 and various fixed values of P 2 -C Q 2 with the GRS densities |1(J . 
The curves refer from top to bottom to P 2 = 0, 0.2 and 1 GeV 2 , respectively. The 
results for the real photon (P 2 = 0) are very similar to the ones in fig. 2 where the 
GRV 7 curves practically coincide with GRS. 

Fig. 7b Our LO distributions as in fig. 7a compared with the ones of SaS ||11| . 

Fig. 8 LO and NLO predictions for the up-quark and gluon distributions of a virtual 
photon 7(-P 2 ) at Q 2 = 10 GeV 2 . For comparison the results for the real photon 
(P 2 = 0) are shown as well. The NLO 'hadronic' contribution in (25) is also shown 
separately. The GRS expectations are taken from ref. [I(J. The DIS 7 results for 
^7(-P 2 ) can fog ea sily converted to the MS scheme with the help of eq. (4), whereas 
g 7 ( p2 ) remains the same in both schemes. The results have been multiplied by the 
numbers indicated in brackets. 

Fig. 9 As in fig. 8 but at Q 2 = 100 GeV 2 . 

Fig. 10 Predictions for the LO effective parton density xf^ p2 \x,Q 2 ) : defined in eq. (29), 
at the scale Q 2 = (Vt*) = 85 GeV 2 and at two fixed values of x. The HI data 
[ [14|j have been extracted from DIS ep dijet production. The solid curves refer to 
our predictions in the theoretically legitimate region p 2 « q 2 = ( P j T et y , whereas 
the dashed curves extend into the kinematic region of larger P 2 approaching Q 2 
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where the concept of parton distributions of virtual photons is not valid anymore 
(see text). 
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